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1. Introduction 

In the original Ehi'enfest scheme for their classification, the order of phase transitions was 
given as that of the lowest derivative of the Helmholtz free energy in which a discontinuity is 
manifest. While first-order solid-liquid-vapour transitions and second-order superconducting tran¬ 
sitions, for example, fit into this scheme, there is only limited evidence for the existence in nature 
of transitions of strictly higher order (i.e., of order three or above). However, on the basis of exper¬ 
iment, recent claims have been made which support the existence of a fourth-order transition in a 
cubic superconductor ||T]] and a theoretical analysis of higher-order transitions was provided in |§]. 

There are many other well known transitions signaled by divergent rather than discontinuous 
behaviour of the appropriate free energy derivatives. By far the most common transitions of this 
type are second-order ones such those associated with symmetry breaking of the Higgs field in 
parficle physics and ferromagnelic fransifions in mefals. The modern classificalion scheme includes 
fhese scenarios and fhe order of a fransifion is now generally faken fo mean fhaf of fhe lowesf 
derivafive of free energy in which any type of singular behaviour occurs. 

Nofwifhsfanding fhe dearfh of experimenfal evidence, higher-order fransifions (especially of 
order fhree) abound in fheorefical models, such as pure spin models |^, spin models coupled fo 
quanfum gravity |Q, and laffice (as well as continuum) foy models of QCD 0. In fhe nofafion 
of spin models, where t = I — T/T). is fhe reduced femperafure {Tc being fhe critical value of fhe 
femperafure variable T) and h = H/k^T is fhe reduced exfernal field (//being fhe field sfrengfh), fhe 
critical poinf is often locafed af (0,0). This may be fhe endpoinf of a line of firsf-order fransifions, 
in which case we assume fhaf t parameferises ifs arc lengfh wifh h in fhe orfhogonal direction. 
We denote fhe free energy by f{t,h) and ifs n* derivatives by fl'"\t,h) and If t or /z 

vanishes, we simply drop if from fhe argumenf. For an m‘*’-order fransifion (t) is continuous for 
n<m — \ while (f) experiences eifher a disonfinuify or is unbounded af t = 0. If is also possible 
fhaf (t) is continuous for n <m' — \ while ^ (t) is singular. The specific heaf, C{t) = (t), 

is fhus continuous if m > 2 and fhe susceptibility, xit) = //, (0> continuous if m' > 2 (a sifuafion 
nof normally possible in a ferromagnef). 

We disfinguish befween fransifions characferised by a disconfinuify in //™^(t) and fhose char¬ 
acterised by divergences. In fhe laffer case, allowing for fhe possibilify fhaf m' is nof necessarily 
fhe same as m, fhe scaling behaviour af fhe fransifion may be described by critical exponenfs: 

( 1 . 1 ) 

In fhe second-order case Qn = m' = 2), fhe exponenfs A and G become fhe usual crifical exponenfs 
a and 7 associated wifh specific heaf and suscepfibilify, respectively. 

Here, some resulfs concerning fhe properties of partition function zeros close fo higher-order 
fransifions are presenfed. The Fisher zeros of sysfems governed by a single fhermal parameter are 
analysed in Sec. ^ where resfricfions on fheir properties are ouflined. The Fee-Yang zeros appropri- 
afe fo fhe field variable are examined in Sec. ^ and a number of scaling relations are derived. These 
recover fhe higher-order relations of 0 in fhe case where m' = m and, when m = 2 , fhey recover fhe 
sfandard Rushbrooke and Griffilhs scaling relations appropriafe fo second-order fransifions where 
fhe specific heaf and suscepfibilify are unbounded. Conclusions are presenfed in Sec. 


244 /2 


PoS(IAT2005)244 



Properties of phase transitions of higher order 


R. Kenna 


Table 1: Impact angles (in the upper half-plane) permitted at a discontinuous phase transition of order m. 
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2. Fisher zeros 


In the situation where the locus of Fisher zeros is linear and can be parameterized near the 
transition point in the upper half-plane by t{r) = rexp(/(/)(r)) - the locus in the lower half-plane 
given by it complex conjugate - the (reduced) free energy and its derivatives are 

f{t)=2Ref ln{t-t{r))g{r)dr , = 2(-l)”^'(« - l)!Re /" dr , (2.1) 

Jo Jo [t-t(r)) 

where /? is a suitable cutoff. 

The difference in the free energy on either side of the transition can be expanded as f{t > 
o)-/(i < 0) = c„i”. For a discontinuous transition, = 0 for n < m, while Cm ^ 0 and the 

discontinuity in the nf-'^ derivative of the free energy is = m\cm- Matching the real parts of 

the free energies across the singular line gives cos n^{r) = 0, so that the impact angle is 


lim(/)(r) = (/> = 

r —^0 


{21 + 1)k 


2m 


for I = 0,... ,m — I . 


( 2 . 2 ) 


The formula ( |2.2D gives a very direct way to read off the order of a transition. A list of permitted 
impact angles in the upper half-plane for discontinuous transitions of various orders is given in 
Table 1. It is clear that while vertical impact is allowed at any discontinuous transition of odd order, 
an impact angle of n/6, for example is only allowed at a transition of order 3 or 9 or 15, etc. On 
the other hand, a discontinuous second-order transition with impact angle 7r/2 is, for example, not 
permitted here. These observations are in accord with various results for first-, second- and third- 
order transitions in [Q], ||^ and Isj] which are associated with impact angles 7r/2 (corresponding 
to / = 0), 7r/4 (/ = 0) and 7r/2 (/ = 1), respectively. The question as to the mechanism by which 
the system selects its /-value is open. 

Assuming analytic behaviour for the density of zeros, it is straightforward to show that the 
necessary and sufficient condition to generate an ni*-order discontinous transition is 


gir) = gor"' ' 


(2.3) 
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where go is a constant. Inserting this into (p?T|), and using ( [2.2[ ), one finds 

A/('”) = 27igo{m— l)!sinm0 = (—l)!27rgo j (2.4) 


which generalizes a well known result in the m = 1 case |^. 

We next consider an ni*-order diverging transition characterised by ( |1.1[ ). To achieve the 
appropriate divergence forit is necessary that 


gir) = gor' 


m—l—A 


Putting this into (|2.1[), one finds fhaf (t) is continuous for n < m, while 




-2go\t\-^r{m-A)r{A)x 


cos (m—A)^ if t < 0 

cos ((m — A)^ +A7r) if ? > 0 . 


(2.5) 

(2.6) 


This general relationship befween critical amplifudes and impacf angles recovers resulfs derived in 
[^] in fhe m = 2 case. If A = 0 if furns ouf fhaf becomes logarifhmically divergenf and 


(?) = 2go(m — 1)!cos (rn^) x 


In I?I if ? < 0 

(In |?| + 71 sin {m(j))) if ? > 0 . 


(2.7) 


When m = 2 fhis recovers known behaviour in fhe second-order case and fhe disconfinuify 
across fhe fransifion is consisfenf wifh ( |2.4{ ). 

A quife nafural origin for finife-size scaling (which is alfernafively and fradifionally viewed 
as a hypofhesis) was offered in by mafching fhe infegrafed densify of zeros on a finile lattice 
fo ifs infinile-volume counferparf. Here, from ( |2.3| ) and (p^), fhis quantify is G{r) r"’ ^ (where 
A = 0 in fhe case of a discontinuous fransifion). For a finife sysfem of linear exfenf L, fhe infegrafed 
densify af fhe y* zero is Gl(|?;|) = {2j — 1)/2L^ @,111, and equating fhis fo G(|?y|) leads fo fhe 
scaling behaviour 

|?y| ~ L~'^ . (2.8) 


In fhe unbounded case, assuming hyperscaling {m — A = 2 — a = vd), fhis recovers fhe sfandard 
form \tj\ ~ for Fisher zeros. In the discontinuous case, where A = 0, ( ^^ yields 


V = 


m 
1 ’ 


(2.9) 


which is a generalization of the usual formal identification of v with 1/d, applicable at a first- 
order transition. Support for this identification comes from the third-order {m = 3) discontinuous 
transition present in the d = 3 spherical model @] as well as in the Ising model on planar random 
graphs if the Hausdorff dimension is used for d 0. 


3. Lee-Yang zeros 

While the Fisher-zero analysis concerns only even exponents, mixed exponents are encoun¬ 
tered in the Lee-Yang case, where both field and fhermal variables enter. Scaling relations between 
such exponents were derived in [||] in the theoretical case where m' = m. These are 

[m — \)A + + G = m{m — , G = p {{m—1)5 — 1) . (3.1) 
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If m = 2, these become the standard Rushbrooke and Griffiths scaling laws, respectively: 

a + 2j3 + 7 = 2, 7 = j3(5-l). (3.2) 


Now, the free energy in the presence of an external field is 


/(t,/i)=2Re / \n{h — h{rf))g{rf)dr, 

JrYL(t) 


(3.3) 


where rYi^{t) is the Yang-Lee edge and the locus of zeros is h{r,t) = rexp (/0(r,t)). The field 
derivative is 


rdx 


(3.4) 


(/r/rvL - 

having used r = xrYL( 0 - We assume that rYL (0 is small near t = 0 so that R/rYL (0 —> °° and, 
when n = m' and h = 0, compare with the limiting scaling behaviour in (| 1 . 1 D to find 


for some function <I> []T3[]. Using n = 1 in (^^ gives, for the magnetization. 


fl^\t,h)=2t-^ rYi.it) 


in'—\ 


Re 


^(x) 


rvL 


liy 


— xe^' 


-dx = t rYhit) 


m! — 1 




^yl(0 


(3.5) 


(3.6) 


Comparison with ( [lI| ) now gives 'R(/r/rYL(0) ~ (/^/'"YL(t)) Eq. (^!^ must become t-independ- 
ent for small t and, under these circumstances, the behaviour of the Yang-Lee edge is 


GS 


(3.7) 


Xyl(0 ~ . 

When m' = 2 and G = 7 , this recovers second-order behaviour. Inserting (^7^ into (^), and 
comparing the /r —> 0 limit of the resulting expression to ( | 1 . 1 [ ), yields the scaling relation 

G 


i 8 = 


{in' — 1)5 — 1 


(3.8) 


h 


V''YL(0y ’ 


(3.9) 


When m' = m, this recovers the Griffiths-type scaling relation ( |3.![ ), derived in [^. 

From ( p.5[ ) and ( p.7[ ), the integrated density of zeros is G{rf) = (r/rYL(f)) 

in which F(x) = /y<I>(x')r/x', enabling (^) to be written as 

for sufficiently small t and where (w) = 2Re F[x )/— x)dx. Its m* temperature deriva¬ 
tive at zero-field is therefore (t) fG(5+l)/({m'-l)5- 

"h an expression which, compared with 

(|lI|), yields the scaling relation 

5 + 1 


A = m-G 


[m' —1)5—1 


(3.10) 


In the case in' = m, (3.8) and (3.1C) recover all four scaling relations derived in [Q], and, in the 
second-order case (m = 2), they recover the standard Rushbrooke and Griffiths scaling laws (3.2). 


244 /5 


PoS(LAT2005)244 





















Properties of phase transitions of higher order 


R. Kenna 


4. Conclusions 

The self-consistency of scaling at higher-order transitions has been examined using partition 
function zeros. While strong constraints on the linear locii (impact angles) of Fisher zeros associ¬ 
ated with discontinuous m^'^-order transitions have been established, the mechanism whereby the 
system selects from m possible locii is an open question and merits further study. For a transition 
characterised by an unbounded m**’ moment, the impact angle is intimately linked with amplitude 
ratios. In both cases, analysis of the impact angle may provide a stable approach to the analysis of 
phase transition order and strength. A complementary analysis of Lee-Yang zeros leads to scaling 
relations between even and odd critical exponents at higher-order transitions of a general type. Fi¬ 
nally, a finite-size scaling analysis shows that the conventional formal identification of V with \/d 
that holds at first-order transitions extends tov = m/d for discontinuous transitions of m* order. 
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